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It has been known that a B = 2 skyrinion is axially symmetric. We consider the Skyrme model 
coupled to gravity and obtain static axially symmetric regular and black hole solutions numerically. 
Computing the energy density of the skyrmion, we discuss the effect of gravity to the energy density 
and baryon density of the skyrmion. 
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' ■ The Skyrme model is an effective theory of QCD constructed by pion fields alone. Baryons are incorporated as 
, topoloeical solitons in this model, called skyrmions, and the topological charge corresponds to the baryon number 
B 0, Ql- Skyrme introduced a hedgehog ansatz for the pion fields to obtain B — 1 spherically symmetric soliton 
^] ' solutions. The detailed analysis for the property of the B = 1 skyrmion as a nucleon was performed in Ref. 
upon quantisation of the collective coordinate. Later Braaten and Carson obtained B = 2 skyrmions numerically and 
showed that B = 2 skyrmions are axially symmetric Interestingly, multi-skyrmions with B > 2 exhibit various 
, discrete symmetries analogously to multi-BPS monopoles 

• The Einstein-Skyrme system in which Skyrme fields couple with gravity was first considered by Luckock and Moss 
, where the Schwarzschild black hole with Skyrme hair was obtained numerically. This is a counter example of the 

• no-hair conjecture for black holes 0- They observed that the presence of the horizon in the core of skyrmion unwinds 
the skyrmion, leaving fractional baryon charge outside the horizon. The full Einstein-Skyrme system was solved later 

O to obtain spherically symmetric black holes with Skyrme hair Q and regular gravitating skyrmions . 

In this paper we shall study regular and black hole solutions of i? = 2 skyrmions with axisymmetry. The recent study 
?-H . by Hartmann, Kleihaus and Kunz showed the Einstein- Yang-Mills-Higgs theory possess axially symmetric monopole 
and black hole solutions Therefore, we shall follow their numerical technique to solve the Einstein-Skyrme 

^ , model with axial symmetry. Computing the energy density and charge density, we shall show how gravity affects to 
the skyrmion. 



II. THE EINSTEIN-SKYRME SYSTEM 

After a long history, the Skyrme model was proved to be an effective theory of QCD by E. Witten in Ref. 0|- At 
low energy the SU{2) x SU{2) chiral symmetry of the strong interaction is broken spontaneously and the Skyrme 
Lagrangian respects this symmetry. The Skyrme model coupled to gravity is defined by the Lagrangian, 

C^Cs + Cg 

where 

Cs = ^.g^''tr([/-ia^[/[/-ia,;7) 

+ ^5^ V^tr( [(7-19^ C/,C/-i9.C/] [U-^dpU,U-^d,U]) (1) 
and U is an SU{2) valued scalar field of pions. f-^ ^ 186 MeV is the pion decay constant. 
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To solve the field equations, ansatz for the metric and Skyrme fields are required. Let us introduce an ansatz for 
the metric given in Ref. ^3 



fdt'+'^idr'' +r^de^) + -y sin^ 



(3) 
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where / — f{r, 6), m = ■m{r, 9), and / = l{r, 9). 

The axially symmetric Skyrme fields can be parameterised by 

U — cos F(r, 9) + if ■ nji sin i^(r, 

with nji — (sin cos niy9, sin sinn(/3, cos 0) and O — 0{r,9). The integer n corresponds to the winding number of 
the Skyrme fields and for S = 2 we have n — 2. 

In terms of F and 8, the Lagrangian takes the form 
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The baryon current in curved spacetime is obtained by taking the spacetime covariant derivative V^, 

247r 

The baryon number then is given by integrating 6" over the hypersurface t = Q, 

B = J drd9dipy^b° 



(6) 



i J drdO {drFdgQ - dgFdrQ) sin 8(1 - cos2F) 
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where {Fq, Qq) and {Fi, 8i) are the values at the inner and outer boundary respectively. In flat spacetime we have 

(i^o,eo) = (7r,0) and (^^i, 8i) = (0, 

which gives B ^ 2. In the presence of a black hole, the integration should be performed from the horizon to infinity, 
which change the values of Fq and allow the B to take a fractional value of less than two. This situation can be 
interpreted as the black hole absorbing a skyrmion. 

The energy density is given by the time-time component of the stress-energy tensor 
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III. BOUNDARY CONDITIONS 

Let us consider the boundary conditions for the chiral fields and metric functions. 

A. Regular Solutions 

At the origin r = 0, the metric functions arc regular, which requires the boundary conditions 

drf{0, 6) = drm{Q, 9) = drl{0, 9) = 0. (9) 

At infinity r = oo, the metric is asymptotically flat. Thus wc have 

/(oo,6i) = to(oo,6I) = /(oo,6i) = 1. (10) 

The boundary conditions for the Skyrme fields at r = 0, cxd arc given by 

F{Q,9) = TT, drQ{Q,9) = 0, (11) 
F{qo, 9) = 0, a,.e(oo, 9) = 0. (12) 

The axially symmetric condition leads to the bomidary conditions on the axes = 0,7r/2 as 

dgfir, 0) = dem{r, 0) = dgl{r, 0) = 0, (13) 

dofi:r, f ) = dom{r, |) = dol{r, \) = 0, (14) 

Sef (r,0) = 0, e(r,0) = 0, (15) 

5,f(r,|) = 0, e(r,|) = |. (16) 

B. Black Hole Solutions 

At the horizon r = r/j, the time-time component of the metric satisfies 

5« = -/(r-/.,^)=0. (17) 

Regularity of the metric at the horizon requires 

m{Th.Q)=l(rh.Q) = ^. (18) 

The boimdary conditions for F{r,9) and 0{r,9) at the horizon arc obtained by expanding them at the horizon and 
inserting into the field equations which arc derived from SCs/6F = and SCs/^^ = respectively. Thus we have 

drF{n„9)=dre{rh,9)=0. (19) 

The condition that the spacetime is asymptotically flat requires 

/(oo, 9) = to(.3c, 9) = ;(oo, 9) = 1. (20) 

The boundary conditions for F and Q at inflnity remain the same as in flat spacetime 

F(oo,6') = 0, 5^e(cx),6») = 0. (21) 

For the solution to be axially symmetric, we have 

dgf{r, 0) = demir, 0) = Oglir, 0) = 0, (22) 

def (r, I) = dem (r, ^) = del (r, ^) = 0. (23) 

Regularity on the axis and axisymmetry impose the boundary conditions on F and as 

deF{r,0) = deF(r,^) =0, (24) 
e(r,0) = 0, e(r,|) =|. (25) 

Under these boundary conditions, we shall solve the Einstein equations and the matter field equations numerically. 
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IV. NUMERICAL RESULTS 

Let us introduce dimensionless coordinate and coupling constant 

X ^ af^r, a = TTGjl. 

Then the free parameters are the horizon xu and the coupHng constant a for black hole solutions while for regular 
solutions the a is the only free parameter. We shall take a = as decoupling of gravity from the matter, effectively 
G = 0. 

In figs, n [3 are the energy densities per unit r and {e — —^fg^TQ) of the regular solutions with a = 0.0, 
3.0 respectively. As a becomes larger, the energy density gets concentrated in a smaller region. Thus the gravity 
influences for the matter to be more compact and denser. The shape remains torus irrespective of the value of a. We 
show the baryon density per unit r and 6 {b = \/'g^b^) in fig. El The shape is similar to the energy density, but its 
dependence on the a is rather small. Note that the energy density and the baryon density both vanish at p = 0. This 
can be readily checked by inserting the Skyrme functions expanded around = into e and b respectively. We found 
that there exists no regular solution for a > 4.79. 

In figs.QlElare the energy densities of the black hole solutions with a = 0.0, 1.5 respectively. As a becomes larger, 
the energy density becomes smaller and sparse. This can be interpreted that the black hole absorbs more skyrmions 
for a larger coupling constant. The shape is slightly distorted in the background of the black hole so that one can see 
the spherically symmetric horizon in the center of the skyrmion. Fig. [HI is the baryon density around the black hole. 
As in the case of the regular solution, the dependence of the baryon density on the value of the coupling constant 
is small. It can be checked that the energy and baryon density vanish at p = in the same manner as the regular 
solution. Inserting the metric functions as well as the Skyrme functions expanded around the horizon instead, one 
can also see that the energy and baryon density vanish at the horizon. 

The domain of existence of the black hole solution is shown in fig. {7\ There exist minimum and maximum value of 
Xh and a beyond which no black hole solutions exist. Therefore the regular skyrmion solutions can not be recovered 
from the black hole solutions by taking the limit of Xh ^ unlike the case of B = 1 ^J. In fig.|Hlis the dependence of 
the baryon number on Xh- One can see that the baryon number decreases as the black hole grows in size. This figure 
confirms that the baryon number is no longer conserved due to the black hole absorbing the skyrmion. 

V. CONCLUSIONS 

We have obtained static axially symmetric regular and black hole skyrmions numerically. As is shown in figs. ^ 
El EI El the energy density is toroidal in shape for the regular and black hole solutions. However, in the black hole 
case, the shape is distorted around the spherically symmetric horizon. The energy density depends on the coupling 
constant. For the regular solution, as the coupling constant becomes larger, the energy density distribute in a smaller 
region and become denser. On the other hand, for the black hole solution, the energy density becomes sparse due to 
the black hole absorbing more skyrmions. The baryon density shows a similar shape with the energy density and its 
dependence on the value of the coupling constant is rather small. However, one can see from fig. El that it depends 
on the horizon size such that the baryon density decreases in increase of Xh- Since there are minimum values of x^, 
the black hole solution can not recover a regular solution by taking the limit of r^, ^ unlike the B = 1 spherically 
symmetric case. 

Obviously it is important to study the stability of our solutions. We expect that the stability analysis may be 
performed by applying the catastrophe theory for black holes with non-linear hair p^ . 

As our further works, it will be interesting to consider skyrmion black holes with B > 3 which have discrete 
symmetries. The inclusion of gauge fields may also be possible to study the interaction between a monopole black 
hole and a deuteron |l3j |. Extending the model to higher dimensions will be also an exciting problem, leading to the 
study of deuteron black holes which may be observed in LHC in future 0, 0, 0| . 
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FIG. 1: Regular solution: The energy density e in cylindrical coordinates p and z with a = 0.0. 
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FIG. 2: Regular solution: The energy density e in cylindrical coordinates p and z with a = 3.0. 
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FIG. 3: Regular solution: The baryon density h in cylindrical coordinates p and z with a = 0.0. 
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FIG. 4: Black hole: The energy density e in cylindrical coordinates p and z with xy^ — 1.0, a = 0.0. 
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FIG. 6: Black hole: The baryon density b with Xh = l-0,a = 0.0. 



